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1 Introduction 

Insulators and superconductors are central concepts in strongly correlated electron materials 
giving rise to high Tc superconductivity, [1]. The less understood region is the region in 
which the insulator-metal transition happens, [2, 3], characterized by novel features of the 
electronic conductivity, [4]. 

The mechanisms that create insulators can be broadly classihed as follows: 
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1. Band gap insulators, where the conduction band is empty, and there is therefore a 
gap that prevents current transport. 

2. Anderson localization, [5], where strong disorder inhibits conduction. 

3. Mott localization, where strong on-site interactions localize electrons, [6]. This has 
been argued to be at work in insulating anti-ferromagnets, [7] as was first pointed out 
by Mott and van Vleck. 

Lately a new further mechanism was argued to lead to insulating behavior: 

4 Momentum-dissipating interactions become relevant (strong) in the IR, and they in¬ 
hibit conduction, [8]. The model used as an example was a holographic model with 
a saddle-point of helical symmetry. The mechanism was generalized to Einstein- 
Maxwell-Dilaton (EMD) theories in [9] and many saddle points corresponding to 
insulators, bad metals and conventional metals were found. 

Holography has provided a new paradigm and a new arena for theoretical models that 
address physics at strong coupling, but in a semiclassical setup (that is typically present 
because of a large-limit involved).^ It is a natural framework to describe quantum critical 
systems at zero and finite density. It is convenient to describe out of equilibrium effects 
as well as conductivity, a major observable in condensed matter. A study of holographic 
ground states has indicated that they are very diverse in their properties, [14], [15, 16], 
[17-19], [20-25]. 

Most holographic systems analyzed at finite density are translationally invariant. Ex¬ 
ceptions also exist, using D-brane defects and magnetic vortices, [26-29], but such systems 
have been much more difficult to analyze so far, although exceptions exist, [30]. The stan¬ 
dard symmetry argument indicates that the real part of the AC conductivity will have a 
S(uj) contribution as in a translationally invariant system a constant electric field generates 
an infinite current. The 5-function can be argued to relate via causality to a ^ pole in 
the imaginary part of the conductivity, via standard dispersion relations. In the context of 
holographic correlators this was analyzed recently in [14]. This 5-function is distinct from 
the one appearing in superfluid/superconducting phases. 

The interaction with momentum dissipation agents has been discussed in rather general 
terms in [31-33]. When the interaction with dissipators is IR irrelevant with respect to 
RG fixed points, a perturbative IR calculation can determine the scaling of the IR DC 
conductivity. When the dissipation is IR relevant, it can change the nature of the saddle 
point, turning the system into an insulator as was first argued in [8]. This is the mechanism 
4 above that leads to insulating behavior. 

There have been several lines of research addressing the breaking of translational invari¬ 
ance in holographic saddle points at finite density and its impact on conductivity. A first 
line of research introduced a holographic lattice imprinted by boundary conditions on the 

^Useful reviews for condensed matter physics applications can be found in [10], [11, 12[ and section 2 of 
[13[. 
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bulk charge contribution, [34-38]. In the regimes accessible to the numerical calculations, 
the lattice perturbation is irrelevant in the IR and it controls to leading order in the IR the 
DC conductivity as predicted in [33] on general principles. 

Another line of approach, [39-44], was to assume an effective action treatment for 
momentum dissipation associated to the breaking of translational invariance and this was 
provided by a massive gravity action. 

More recently a third approach to breaking translational invariance and therefore in¬ 
troducing momentum dissipation was introduced in [45] and [46]. In [45] the phase of UV 
relevant, charge neutral, complex scalars was used in order to break translational symmetry 
on the boundary. The construction of [46] used bulk fields which are perturbatively without 
a potential on AdS. In order to relax momentum, the massless scalars were given a linear 
dependence on the spatial coordinates of the boundary. 

The formula obtained for the DC conductivity in all the studies above is a sum of two 
contributions: 

The first term, has been interpreted, [47], as a pair creation contribution as it persists 
at zero charge density. For the RN black hole it is a constant proportional to the inverse 
of the bulk gauge coupling constant that counts the relative density of charge-carrying 
degrees of freedom to the neutral ones in the strongly-coupled plasma. More recently, it 
was verified in [48] that the first term in (1.1) is the electric conductivity in the absence 
of a heat current. The interpretation as a pair creation term is then natural since charged 
pairs are created with zero total momentum and therefore not contributing to a net matter 
flow. It is fair to say however that this pair production interpretation, although appealing is 
not in agreement with other related expectations from weak coupling, like the fact that we 
expect it to be exponentially suppressed in the chemical potential as this gaps free fermions. 
At finite density, this term persists and can be powerlike in the temperature in scaling IR 
geometries. We will henceforth call this term the charge conjugation symmetric (ccs) term 
following J. Zaanen’s suggestion. 

The second contribution in (1.1) is due to the effects of dissipating momentum. When 
translation-breaking operators are irrelevant, the system is expected to be metallic and this 
term is expected to give the leading contribution to the DC conductivity. Then, a descrip¬ 
tion of momentum relaxation in terms of the memory matrix formalism is appropriate and 
shows that the conductivity takes a Drude-like form, though no quasi-particle description 
is assumed [33]. 

This general form of the DC conductivity was seen already in pure metric backgrounds 
in [47] and was generalized to dilatonic backgrounds in [17]. In both cases, as the gauge 
field action is the DBI action, (1.1) is replaced by 

<yDC = V^(u|f^)2 + (ag^*)2 (1.2) 

giving results compatible with (1.1) in the regimes where pair creation or diffusion dominates 
the conductivity. In general, we expect a nonlinear formula that reflects the bulk action of 
the gauge field. In the probe DBI cases the momentum dissipation is due to the fact that 


- 3 - 



charge degrees of freedom are subleading compared to uncharged ones. This means that 
there is a momentum conserving (5-function but its coefficient is hierarchically suppressed 



In [17] it was observed, based on (1.2), that for running scalars other than the dilaton 
and in 2+1 boundary dimensions, the drag DC resistivity, when it dominates, is proportional 
to the electronic entropy. This is a general property of strange metals where both the 
measured electronic entropy and resistivity are linear in temperature. This was extended in 
[49] to more general cases using the massive graviton theory, and most importantly provided 
a kinetic explanation for the correlation suggesting a more general validity. 

The general properties of holographic conductivity were further corroborated recently 
by a careful study of the holographic current-current correlators, and the associated prop¬ 
erties of their poles in the complex plane, [50]. 

In most holographic cases, the translationally invariant systems are good conductors 
and we will call them “holographic metals”. There have been also systems that have a 
discrete spectrum in the current-current correlator and they are therefore candidates for 
insulators. At zero density such systems were described in [51-54]. 

More interestingly, at finite density a large class of systems were found in [17, 18] 
that had a discrete spectrum for the current-current correlator (as well as the stress-tensor 
correlators). One particular case in a different supergravity was independently observed in 
[55] and a new one (the two-charge STU solution in AdS^) is presented in this paper. Such 
holographic saddle points have very peculiar properties. In particular, the spectra of the 
stress tensor and current-current correlators were discrete. The current-current correlator 
had a (5-function at zero frequency because of translational invariance. Finally, the IR 
endpoint of the bulk geometry does not affect IR physics but UV physics. Moreover, 
although the analysis of [17] was based on “bottom-up” theories, the example of [55] and 
the one studied here are found in M-theory compactifications and are both therefore “top- 
down". 

Interestingly, the thermodynamics for such systems is reminiscent of that of Yang-Mills 
in four dimensions, [56, 57]. At finite temperature, up to a transition temperature Tc, the 
dynamics is temperature-independent to leading order in the large-A^ limit. However there 
is a first-order transition at Tc to a new phase that is strongly coupled, liquid-like and 
gapless. However, in these theories, the saddle points are translationally invariant and the 
DC conductivity is strictly infinite because of the (5-function. 

The purpose of this work will be to introduce momentum dissipation in the holographic 
theories in question and study the ensuing holographic saddle points. We will model momen¬ 
tum dissipation by using massless scalars with translationally non-invariant backgrounds, 
[45, 46] and choose it to be irrelevant in the IR. Our results are as follows: 

• In EMD theories, as in (1.3) with no U(l) symmetry breaking (!+((/>) = 0) and with 
momentum dissipation we analyze the nature of the zero-temperature current-current 
correlator for theories that without momentum dissipation had a discrete spectrum. 
We find that the zero-temperature spectrum remains discrete but as expected the zero 
frequency (5-function disappears. Such saddle points are insulators with a hard gap 


- 4 - 


and a discrete spectrum. They share similarities with insulators of type 1 above, with 
the difference that there, electrons are weakly-coupled and the underlying spectrum 
of possible states is not-only gapped but has bands. Therefore the only similarity 
with type 1 insulators is the presence of a gap. In our case this is driven by strong 
interactions in the underlying theory, and in that respect it resembles type 3 (Mott) 
insulators although the detailed mechanism seems different. 

Such ground states share many similarities with nuclear matter at finite charge density, 
[58]. 

We believe that this mechanism of insulation is new and it should be number five on 
the previous list. To substantiate this claim however further work is needed and we 
comment on this in the outlook section. 

• To corroborate the previous claims a subtle and detailed analysis of the T = 0 current- 
current correlator, and the conductivity are required. The reason is that at T = 
0 a zero frequency 5-functions can re-emerge as it happens in metals, because the 
momentum dissipating effects are IR irrelevant and therefore shut-off at T = 0. Our 
analysis indicates that in the saddle-points in question this does not happen. It is 
interestingly correlated with another feature of non-trivial quantum critical saddle 
points: the presence of the mild IR (resolvable) singularity and its resolution. As was 
first argued in detail in [17], sometimes the calculation of a correlator in the presence 
of a resolvable singularity is well defined and the result independent of its resolution. 
We called such cases, “repulsive" or “holographically well-defined". We show here that 
there is a clear correlation between such holographically well defined cases and the 
absence of the zero frequency 5-function. 

• We construct examples of such saddle points numerically using a concrete bulk grav¬ 
itational theory and a complete flow from the UV to the IR. As the zero temperature 
geometry is numerically challenging we approach it from the unstable black hole 
branch. Our findings corroborate the asymptotic analytical tools used to analyze the 
zero temperature saddle points. 

• We further consider effective holographic actions with U(l) symmetry breaking (spon¬ 
taneous or explicit) as in (1.3) with (IT((/>) ^ 0). We consider fractionalized IR fixed 
points of the type found in [19], with a discrete spectrum for the current-current cor¬ 
relator. We add momentum dissipation as in the previous case and verify analytically 
that the nature of the spectrum of the current-current correlator remains similar: it 
is discrete. In the case where the U(l) symmetry breaking is spontaneous, there is 
always a 5-function at zero frequency due to the superfluid pole. Such saddle points 
describe superfluids in the presence of the breaking of translation invariance, with a 
discrete spectrum and a hard gap above the superfluid pole. Their behavior in this 
respect resembles supersolids, ^ although in the supersolids discussed in the literature 

^Supersolids (see [59] for a review) have been anticipated theoretically, [60] and have been described 
in several contexts, [61-67]. Experimental evidence for supersolids has been claimed in [68] . Realizations 
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the translation invariance breaking is expected to be spontaneous. We may call them 
“charge supersolids” although the spectral density of the spin-2 (shear) part and the 
spin-0 (bulk) part of the stress tensor is also discrete and gapped in these saddle 
points (as detailed in section 7). 

The rest of this paper is organized as follows: 

• In the next subsection the main gravitational (bulk) Lagrangians are introduced and 
the procedure for constructing the extremal IR geometries (along the lines of [17-19]) 
is described. 

• In section 2, the equations of motion for the holographic saddle points are analyzed 
and an analytic forms of the asymptotic solutions are obtained and studied. 

• In section 3, the conductivity is studied analytically for scaling solutions of the EMD 
theories introduced in section 1.1. Such solutions contains both the charged and 
neutral hyperscaling violating solutions of [17, 18]. 

• In section 4, the conductivity is analyzed in the presence of momentum dissipation. 
In particular it is shown that the nature of the charged spectra do not change in the 
relevant theories. 

• In section 5, a further analysis of the DC conductivity substantiates the claim that 
such systems have properties reminiscent of strongly coupled insulators and super¬ 
solids. 

• Section 6 contains the numerical results in concrete holographic models that substan¬ 
tiate the analytical conclusions obtained in previous sections. 

• In section 7, we consider the spin-2 and spin-0 parts of the stress-tenror two-point 
function. 

• In section 8, we present our conclusions and outlook for the present work. 

• In appendix A, we analyze in detail the equations of motion. 

• In appendix B, the potential for the scalar, used for the numerical evaluations is 
analyzed in detail. 

• Appendix C contains an overview of the properties of charged black holes known from 
supergravity, and compares them with the systems described in this paper. 

• Finally in appendix D, we analyze zero density systems that are gapped, by applying 
the same method. 

in cold atoms have been also pnt forward, [69]. Most recent experiments however cast doubt on the first 
experimental observation, [70] although there are further experimental claims for verifications. 
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IR fixed points 



Figure 1. The possible IR geometries from systems with a U(l) symmetry which may be broken 
or unbroken. The running scalar cases correspond to those that may have a gapped geometry. 


1.1 Holographic finite density systems at quantum criticality 

We will consider a general Einstein-Maxwell-scalar action in fonr dimensions with a massive 
gange field (mEMD) 


5 = y 






2=1 


(1.3) 


where (t = 1, 2) are massless scalars responsible for the breaking the translational 
invariance. We will eventnally reqnire the last term to be irrelevant in the IR. Moreover 
these scalars will take a simple form, namely tpi = kiXi (no snm on i) in order to break 
translational invariance. 

The scalar (p is dnal to the most important scalar relevant operator of the system 
that drives a non-trivial RG flow from the UV to the IR. We assnme that Z{p), 

and Y{(p) have exponential asymptotics as in snpergravity systems. To obtain an 
asymptotically AdS geometry, V{p) mnst have a finite extremnm. The detailed analysis of 
the IR geometries from (1.3) withont the axions was given in [19]. 

The possible IR scaling geometries are snmmarized in fignre 1 where we nse the ter¬ 
minology of reference [19] to characterize the holographic critical scaling geometries. In 
particnlar, the two sides of the diagram refer to the bnlk action nsed. On the right-hand 
side, the bnlk action has the U(l) symmetry nnbroken (or the dnal gange boson is mass¬ 
less). On the left-hand side the symmetry is broken (either spontaneonsly or explicitly) and 
the gange boson is massive. As we stndy a finite density system, the gange field is always 
non-trivial in the bnlk backgronnd solntion. There are two possibilities: 


• The IR-neutral solution is a solntion where the gange field vanishes fast enongh and 
therefore does not affect the backgronnd solntion in the IR. This does not necessarily 
imply that the electric finx vanishes in the IR. 

• On the other hand in IR-charged solutions the gange field is affecting non-trivially 
the solntion in the IR. 
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Figure 2. A typical relation between the temperature and (j)fi = (j){Th) for the gapped geometry. 
There is a minimal temperature, and the ground system can be obtain by taking the extremal limit 
T —>■ oo. The dashed line shows the relation between T and 4>h when the axions are added, while 
the solid line shows the case without axions. The parameters are fc = 1/2 and A = 1/2. 

The scalar cf) can either be a constant or rnnning in the IR. We are interested in the 
case when (p is rnnning, which leads to the IR geometry 

(1.4) 

where z is the Lifshitz scaling exponent, and 9 is the hyperscaling violation exponent. We 
will change f to r later. The exponents {z, 6) do not nniqnely characterize the solntion. 
Beyond the metric the behavior of the gange held is important to characterize the solntion 
and there is another critical exponent (p that was introdnced in [19] and characterized in 
more detail in [71].^ 

In the case where the U(l) symmetry is broken and therefore there is a nontrivial mass 
term iT((/>)A^ for the gange held in the bnlk, we need to nse two more concepts. 

• If the symmetry breaking term does not contribnte to the leading order of the IR 
geometry, then the electric hnx in the IR 

j *[Z{(P)F] ~ (1.5) 

is constant, and the system is in a fractionalized phase. 

• If the symmetry breaking term contribntes to leading order to the IR geometry, then 
the electric hnx in the IR is zero, and the system is in a cohesive phase. 

We will mainly stndy the fractionalized phase in this work. See [74, 75] for details 
abont the fractionalized and cohesive phases. 

We now focns on the IR charged, fractionalized phase, and briehy mention the other 
two hyperscaling violating cases, which are cohesive phases. In the fractionalized phase, 

somewhat similar exponent was also introduced for different reasons in [72] and also in [73]. 
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(b) 






Figure 3. Schematic plots of the effective potential and the conductivity with a gap. (a) The 
spectrum is discrete, (b) The spectrum is continuous with a gap. 

the metric exponents (z, 9) can be obtained in terms of two parameters in the Lagrangian 
( 7 , 6). These are dehned in terms of the asymptotic form of the potential functions in (1.3) 
in the far IR of the solution as 

V{cl)) ~ , Z{(t)) ~ , W{(t)) ~ e"'^. (1.6) 

The Gubser criterion for the singularity that is present in these backgrounds is implying 
[17, 76]^ 


37^ - 27(5 - (5^ + 4 > 0, 7 ^ - 7(5 + 2 > 0, (5^ - 7(5 - 2 < 0. (1.8) 


When these inequalities are satished, the naked singularity is expected to be resolvable. 
Moreover, we require that the singularity is “repulsive". In practice this means that the 
Sturm-Liouville problem (associated to the conductivity calculation) be well dehned, i.e., 
there is only one normalizable solution in the IR limit, among the two linearly independent 
solutions. Under the above conditions, the naked singularity is expected to be resolvable^ 
and moreover the conductivity does not depend on how this resolution is done. 

^This is stronger than the null energy condition: (2 — 6){2z — 2 — 0) > 0, {z — 1){2 + z — 9) > Q [77]. The 
Gubser criterion in terms of z and 6 is [18] 


z + 2-e 
2Z-2-9 


> 0 , 


2-61 + 1 
2Z-2-9 


> 0 , 


2-1 

2Z-2-9 


> 0 . 


^Either by embedding in a higher-dimensional geometry or by a stringy resolution. 


- 9 - 


(1.7) 





















The conductivity for holographic models satishes a Schrodinger-like equation. If the 
Schrodinger potential V only allows bound states, then the spectrum is discrete. This 
implies that the poles of the conductivity have a discrete pattern and the conductivity is 
given by a discrete sum over the poles. See hgures 2 and 3. 

The two-dimensional parameter space ( 7 , 5) will be characterised on the basis of con¬ 
ductivity. We plot the parameter space for ( 7 , S) when 1/ —)■ 00 in the IR in hgure 4 below. 
If y —)■ 00 in the UV also, then Re[fT(a;)] has discrete spectrum only; if R is a positive 
constant in the UV, then Re[cr(t<;)] has a hard gap, may have a few low-lying poles and 
continuous spectrum above this. The parameter space when the conductivity has a gap 
will be identihed, and the corresponding geometry at hnite temperature will be numerically 
constructed. 

The DC conductivity for a translationally invariant system at hnite density has a h- 
function at cu = 0. This is because in a translationally invariant system charge at rest can 
be boosted to arbitrary velocities implying inhnite DC conductivity. 

Momentum dissipation can resolve the (5-function. In this work we will choose as 
the mechanism for momentum dissipation a distribution of string charge, described by an 
inhomogeneous solution of scalars without potentials,® [45, 46]. The huctuation equations 
that determine the conductivity involve two coupled equations that couple the gauge held 
huctuation to the axion huctuation responsible for the momentum dissipation. When the 
effects of momentum dissipation on the background solution vanish in the IR, ^ we show that 
the qualitative behavior of the conductivity is the same as without momentum dissipation, 
except that now the (5-function peak has disappeared. The gap in the conductivity, and the 
discrete spectrum of poles persists. 

To determine whether the system is a metal or insulator, we must calculate the DC 
conductivity. The DC conductivity of holographic systems has two types of contributions. 
One exists and is hnite without the momentum dissipation, and is also non-trivial at zero 
charge density. It has been interpreted as due to charge pair-production albeit in a critical 
massless system, [47] and this interpretation was reinforced in [17] and [48], although this 
interpretation is still tentative. 

The other term originates from the momentum dissipation. If there are several sources 
of momentum dissipation then there are several such contributions that add-up. As seen 
in special cases in [9], this second term has the structure , where Q is the charge 

density and Tj are the various diffusion rates. 

In the class of theories that we study the dissipative term always dominates. This 
is true for both gapless and gapped geometries. In each case the system can be either a 
conductor or an insulator. We will later obtain the parameter space for which the system 
is a conductor or insulator. 

We hnd theories that in the presence of momentum dissipation have a current two- 
point function whose imaginary part is a discrete sum of poles on the real axis. This gives 
a conductivity whose real part is a sum of (5-functions. This behavior persists at all non- 

®These are known as axions in string theory. 

^This is tantamount to the dissipation dynamics being IR irrelevant. 
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zero temperatures T < Tc up to a fist order phase transition in which the system jumps 
into a “plasma" phase described by a hairy black hole. The current correlator at leading 
order in 1/Nc is independent of temperature in the low temperature phase 0 < T < Tc- 
The system in this phase transport-wise resembles a Mott insulator.® This behavior was 
conjectured hrst in [17] after translationally-invariant solutions with a discrete spectrum in 
the current-current channel was found.® 

Alternatively, we analyze holographic systems that are effective theories with a broken 
U(l) symmetry in the presence of momentum dissipation (arising from broken translational 
invariance), along the lines of [19]. We again use axion solutions that break translation 
invariance. We hnd holographic superfluid holographic ground states, whose two-point 
function of the current has a discrete set of poles including the pole at zero frequency 
guaranteed by Goldstone’s theorem. These systems are superconductors with a discrete 
spectrum. In many respects they resemble supersolids, [59, 69]: Translational invariance 
is broken, there is a pole in the conductivity at zero frequency, and they have a discrete 
spectrum of lattice vibrations. A main difference is that supersolids are crystalline while 
the holographic system is (not-surprisingly) continuous. Like the Mott-insulator-like ground 
states described above, here also there is a low temperature supersolid phase up to a hnite 
transitional temperature, T^- To leading order in 1/W the spectrum and conductivity are 
temperature independent. The system enters a liquid plasma phase (black hole) above Tc- 

To determine whether such IR ground states can be reached by RG flows from the 
UV we solve the hnite temperature equations of motion and then take the near-extremal 
limit. We explicitly construct two systems with hyperscaling violating geometries in the 
IR. The near-extremal scaling behaviors are obtained and matched with the analytic results 
from the near-extremal black hole. This gives rise to complete holographic ground states 
realizing the behaviors (Mott-insulator and supersolid) described above. 


2 Constructing the extremal geometry 


The detailed analysis of the IR geometries from (1.3) is given by [19]. We will construct 
the hyperscaling violating geometries from a Einstein-Maxwell-scalar system with VL = 0. 
This is sufficient for us to discuss the fractionalized phase, in which the symmetry breaking 
term is subleading in the IR. 

We want to obtain the extremal geometry of the following system 


S= 




( 2 . 1 ) 


where 

Z{c^) = !/(</.) = Voe-^^ + + • • • . ( 2 . 2 ) 


We are looking for the extremal geometry when (()—)• oo in the IR. We assume that the 
leading term in V{4>) in the IR is the hrst term and Vq < 0 is required. 

®Mott insulators are defined as insulators due to strong repulsive on-site interactions. In this paper we 
use a more liberal and more general definition: insulators due to strong interactions. 

^Another example of a translationally invariant holographic finite density system with a discrete spec¬ 
trum was found independently in [55]. 
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The metric ansatz is 


= 1 I - 


q(r) , 9 ,9 

—rdt H——- + dx + dy'^ 
h[r) g[r) 

The AdS boundary is at r = 0, and the IR limit is r —>■ 00 . The IR geometry is 


1 




26 


ds^ = ^ “1 —dr + r2-«-h dx^ + dy 


ho 


90 


(2.3) 


(2.4) 


where z is the Lifshitz scaling exponent, and 6 is the hyperscaling violation exponent. 
The equations of motion are 


2r‘^gh4>” — Arghcj)' + 2r‘^hg'(f>' — r‘^gh'(j)' — 2hV'{(f)) + r^h? Af Z' {(j)) = 0, 

{Vhz{^)A[y = 0 , 

Arg' - 12g - 2V{(p) - r'^hZ{(p)A[^ - r'^gcj)''^ = 0, 

2h' - rhcj)'"^ = 0 . 


There is a conserved (extremality) charge 

Q-rh[z{4,)A,A',-y{f) 

which is zero for the extremal geometry [16]. Equation (2.5b) can be solved as 

VhZ{ 4 ))A[ = -p , 


(2.5a) 

(2.5b) 

(2.5c) 

(2.5d) 


( 2 . 6 ) 


(2.7) 


where p is the charge density. We may use this equation to eliminate At in (2.5a), (2.5c), 
and (2.5d), and obtain three coupled equations for (j), g, and 
The near horizon expansion for the functions is as follows 


e-^ = r-(/o + hrP + • • • + fir^ + + •••), 

At = r^“^(ao + dir^ + d2r 


2/3, 


2/3, 


+ Oir" + a2r^" + • 
• + 513-“ + 523-^" + 


g = r 2-6 (^Q + + g 2 r 

h = r ^ 2-s (/jQ _|_ + h2r‘^h - 1_ H-), 


(2.8a) 

(2.8b) 

(2.8c) 

(2.8d) 


where the leading terms are from the hyperscaling violation geometry, the terms involving 
(5 are the modes of the linear perturbations and the terms involving a are from the second 
exponential in the scalar potential. There can be more powers if the potential has more 
exponentials. 

3° The relation between r and f in (1.4) is The advantage of the coordinates (2.3) is that it 

is the Poincare coordinates in the UV, and thus more convenient for numerical calculations. 

^^We can use the conserved charge Q to determine the chemical potential fi. At the horizon, Q = Ts, 
where T is the temperature, and s is the entropy density; at the AdS boundary, Q = —^J-p + 3e/2, where 
e is the energy density determined by (6.6) below. However, if we work in the grand canonical ensemble 
(fixing n), we need to solve the four equations. 
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(2.9) 


The exponents k and are defined in terms of the coordinates (1.4) [17]: 




A 


The relation between k, ( and k, (( is 

2k 


K = 


t ~ r 


C -2 = 




2{C-z) 


( 2 . 10 ) 


2 -e' ^ 2 -e 

After we substitute the above expansion to the equations of motion, we obtain equations 
order by order in the perturbations. At the zeroth order, we can solve /t, C, z, 6, go, and 
oq as follows. The parameter ho is hxed by requiring h = 1 in the UV. Thus, the only free 
parameter in the IR is /q. There are two sets of consistent solutions. 


IR charged solution, [17]: 


z = 


7 ^ + 2^5 — 35^ + 4 

72-52 ’ 


e = 


AS 


K = 


90 = 


al = 


7 — 5’ 


K = 


7 + 5 

/o“'(7-5)Vo 


7 + 5’ 

C = 0-2, 


( 2 . 11 ) 


2(372 — 2 jS — 52 + 4)(72 — 75 + 2) ’ 

fo 


2/o"^"''(7-5)^(52-75-2)Ro 


^0(372 — 275 — 52 + 4)2(72 — ^5 -|- 2) 

The powers fSi arise from perturbing the hyperscaling violating geometry, 
4 + (7 — 5)(37 + 5) 


do = 0 

d± = 


2(7 - 5)2 

4 + (7 — 5)(37 + 5) ± •\/(4 + (7 — 5)(37 + 5))(36 + (7 — 5)(175 + 7(19 + 8(7 — 5)5))) 


2(7-5)2 

In the IR expansion, /3j can be taken as either /3o or 
IR neutral solution: 

252 


( 2 . 12 ) 


z = 1, 0 = 

K = —25, 

90 = 


52 - 1 ’ 

25 


K = 


52 -r 


C = 


2(52 - 75 -1) 

5^^^1 ’ 


(2.13) 


fa‘Vo 


2(3-52)’ 

We will only focus on the IR charged solution in this work. 

The terms involves a are from other exponentials in the potential of the scalar. At 
the hrst order, we can solve /i, ai, gi, hi, and a. For the hrst set of solutions (2.11), the 
solution for a is 

4(5+ 5i) 


a = 


7 — 5 


(2.14) 


we write the numerator of /3± as a ± \/h, we have a > 0 and — b = 8(37^ — 275 — 5^ + 4)(7^ — 75 + 
2)(5^ — 75 — 2) < 0, according to the Gubser criterion; consequently, we always have /3- < 0 and /?+ > 0. 
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For the second set of solutions (2.13), the solution for a. is 

a =-25(5 + (5i). (2.15) 

Note that a depends on the UV completion of the system. The value of a gives a further 
constraint for ^ and 0 by requiring a < 0. 

The asymptotic behavior of the scalar field cf) near the AdS boundary is 

4> = -) + (?ifer^“^(H-), (2.16) 

where (j)a is the expectation of the scalar operator dual to (j), and is its source. When A 
is an integer, there may be log terms. We discuss the choice of the potential V{4>) to avoid 
the log terms in appendix B. 


3 Conductivity with translational invariance 


3.1 The IR charged solution 

We will focus on the conductivity calculated from the extremal geometries. We are in¬ 
terested in the case when the conductivity has a discrete spectrum for the hyperscaling 
violating geometries. To obtain the conductivity, we perturb the system (2.1) around the 
solution to (2.5) by 5Ax = ax{T)e~^‘^^ and 5gtx = ■ Generically, for the metric 

ds'^ = —D{r)dt^ + B{r)dr^ + C{r){dx‘^ + dy^), (3.1) 


the equation for ax after eliminating gtx is, [17] 

/ / 


ID 

b' 


Z\l^a' 1 + ( Zxl^uj^- 


'B 


Z^A?\ 

7^) 


Qx = 0 . 


(3.2) 


After we impose an appropriate boundary condition in the IR (to be discussed below), the 
asymptotic behavior in the UV is 


ax{r) = -h a^r + ■■■ , 
and the conductivity is calculated from^^ 


G = 


a7 


ICO 


(3.3) 


(3.4) 


We do not need to know ax explicitly to decide whether the spectrum of Re(a;) is gapped/gapless 
and/or discrete/continuous. 

After a change of variables by 


dx — (^■^) 

Z — Zq at the AdS boundary, then there will be an extra factor in G so that G — Zq^^. We can 
set Zo ~ 1 hy rescale A^. This amounts to a suitable change of units of measuring the boundary charge 
density. 
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we can obtain a Schrodinger equation 

d O^X -r~r / J.\ 9 

- + V{C)ax = o?ax. 

The potential is given by 

~ _ ZAf DB'Z' D'Z' DZ'"^ DZ" 

~ ~l3 452 Z 4BZ ~ Ib^ 2BZ' 


(3.6) 


(3.7) 


where the prime is denoting a derivative with respect to r. If asymptotes to infinity in 
both the IR and the UV, then the Schrodinger eqnation can only have bonnd states. 

To analyze the behavior of V in the IR, we nse the hyperscaling violating geometry 
(2.4) as backgronnd: 


_ 4(-2 — 1 ) _ n 

D = r 2 -e 


We change the variable r to ^ as 


B = r 2 -i 


;-2 




^ V'^ 2 — 0 _ 2 ^ 

A = 2—11-y.2-9 

ffo 2z 

The asymptotic Schrodinger potential is given by 

z/2-1/4 


where 


V = 




(7 — h )(37 + 5(5) + 12 3z — 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


2[(7 - <^)(7 + 3h) + 4] 2z 
The expression of v in terms of z and 0 holds for both the fractionalized phase and the 
cohesive phase when the cnrrent is dominant in the IR.^'^ 

The bonndary condition in the IR depends on whether the IR limit is at 00 or 
0. This is directly related to whether in the extremal limit the black hole solution has 
T — 7 > 0 or T —)• 00 in the near-extremal geometry. The relation between the entropy and 
temperatnre is 


Comparing (3.9) and (3.12), we can distingnish three different cases as follows: 


(3.12) 


(1) The IR limit is at 00 , which happens when > 0. In this case, the extremal 
limit of the small black hole branch is at T —)• 0. Moreover, the Gnbser criterion implies 
that we have > ^. To leading order in the IR expansion, the solntion for ax with the 
in-falling wave bonndary condition is 

~ (3.13) 


The cnrrent-cnrrent correlator is gapless in this case. 

use two different terms to describe two different effects. We call the current “dominant" in the IR 
if the gauge field is both non-trivial in the IR and contributes to the leading order IR background solution. 
This is distinct from the term “relevant" current that indicates a current perturbation of the IR solution 
that is relevant in the RG sense. 
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(2) The IR limit is at ^ 0, which happens when < 0- In this case, the extremal 

limit of the small black hole branch (that is now thermodynamically unstable) is at 
T —> oo. Moreover, the Gubser criterion implies that v < 0. To leading order in the IR 
expansion, the general solution for dx is 

dx = Civ^J_,(a ;0 + ^^ 2 V^J.(w 0 ~ Ci ^ ^ 1 / 2 +.^ 

The first linearly independent solution in (3.14) is always normalizable: 

dxl'^d^ < 00 . (3.15) 

The second linearly independent solution in (3.14) is normalizable when \u\ < 1, and is 
non-normalizable when \v\ > 1. If there are two normalizable solutions, the boundary 
condition in the IR depends on how the singularity is resolved, and thus the calculation 
of the correlator is unreliable and can only be fixed when the singularity is properly 
resolved. 

When |z^| < 1/2, the current-current correlator is gapless; when \v\ > 1/2, the current- 
current correlator is gapped, if the Schrodinger potential in the UV is also positive, 
which happens for scaling dimension 1/2 < A < 2 .^® 

(3) The IR limit is at a constant non-zero which happens when 2 ^ —t 0. In this case, 
the extremal limit is at a constant T, and the Schrodinger potential is a constant Vq. 
The general solution for ax is 

ax = (3.16) 

When > Vq, the first solution describes the in-falling wave with the oj ^ u + ie 
prescription. When < Vq, both solutions are real and normalizable. As the back¬ 
ground extremal solution is singular, this situation seems to be in the holographically 
non-well defined class we have discussed earlier: we need to impose an extra boundary 
condition at the singularity in order to obtain a unique solution. However, in this case, 
analyticity (in oj) of the correlator makes the solution unique. Indeed, if we analytically 
continue the > Vq solution to < Vq, the solution for ax is unambiguous for all uj. 

When Vq < 0, the current-current correlator is gapless; when Vq > 0, the current-current 
correlator is gapped, if the Schrodinger potential in the UV is also positive.^® 

If the IR is at ^ — 7 - 00 and \v\ > 1/2, the potential V is divergent in the IR. To 
summarize, the conditions for a gapped background charged spectrum are 

• Gubser criterion: (7 — (5) (37 -|- 5) + 4 > 0, (7 — + 2 > 0, (7 — <5)5 -|- 2 > 0 [17]. 

[17] where an analogous condition was first derived it was assumed that near the UV fixed point at 
0 = 0, Z((/)) = 1-1- Here we assume the more general case = 1-1- 

^®It is fair to say that such cases are degenerate and generalize a similar case that was discussed in [56, 57[. 
Further and deeper analysis is necessary in order to find out waht happens in such cases. 
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Figure 4. This figure refers to the IR charged solutions. The regions A, B, and C are the 
parameter space constrained by the Gubser criterion. We assume that the conditions detailed at 
the end of section 3.1 in the UV for a gapped geometry are satisfied. In the left plot the axes are 
given by the gravity action parameters ( 7 ,( 5 ). In the right plot the axes are given by the metric 
critical exponents {z,9) that are related to ( 7 ,5) as in (2.11). In region A (yellow), the extremal 
limit is at T —>■ 0, and the current-current correlator is gapless. In region B (red), the extremal limit 
is at T —> 00 , and the current-current correlator is gapped. In region C (green), the extremal limit 
is at T —> 00 , and the current-current correlator is gapless. Region D (enclosed by blue boundaries) 
is holographically unreliable. 


• Thermodynamical instability of the small black hole branch (^y^ < 0): {j — S){'y + 
3(5) + 4 < 0. 

• The Sturm-Liouville problem is well-defined: |z/| > 1 (stronger than |z/| > 1/2). 

Figure 4 is a summary of the above discussion. The parameter space when the current- 
current correlator can have a gap is region B. 

The Schrodinger potential V near the AdS boundary was analyzed in the Appendix 
of [78] and their conclusions are as follows. The asymptotic behavior of cj) near the AdS 
boundary is (2.16). When both (/)a and cjih are nonzero, V in the UV is divergent when 
A 7 ^ 1 and A 7 ^ 2, and is constant when A = 1 or 2. When the source is zero {<pb = 0), U 
in the UV diverges when 1/2 < A < 1 and 1 < A < 2, is constant when A = 1 or 2 , and 
vanishes when A > 2 (the unitarity bound is A > 1 / 2 ). 

3.2 The IR neutral solution 

When the current is subdominant in the IR we must use the IR neutral solution solution 
(2.13). To leading order in the IR geometry, the effective potential for the Schrodinger-like 
equation is 

V{x) = (3.17) 
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Figure 5. This figure refers to the IR neutral solutions. The regions A, B, and C are the parameter 
space constrained by the Gubser criterion. We assume that the conditions detailed at the end of 
section 3.1 in the UV for a gapped geometry are satisfied. In the left plot the axes are given by 
the gravity action parameters ( 7 ,< 5 ). In the right plot the axes are given by the critical exponents 
of the metric 9 and gauge field C, defined in (2.10) and given in (2.13). In region A (yellow), the 
extremal limit is at T —>■ 0, and the current-current correlator is gapless. In region B (red), the 
extremal limit is at T —>■ oo, and the current-current correlator is gapped. In region C (green), 
the extremal limit is at T —>■ oo, and the current-current correlator is gapless. Region D (enclosed 
by blue boundaries) is holographically unreliable. In region E (enclosed by black boundaries), the 
conductivity has a 5-function at w = 0. 


where 


52 - 275 - 1 _ c - 1 

2(52 - 1) “ 2 


(3.18) 


The analysis is similar to the case with a dominant current, except that the values z, 6 and 
( are given by eq. (2.13). 

When the IR is at ^ 0, the solution for is 


ax = 


(3.19) 


where t'o is defined in (3.18). When \i^q\ < 1, both solutions are normalizable, and this 
is a holographically unreliable case as argued above in the sense that an extra boundary 
condition is necessary at the IR singularity. Therefore we require t'o > 1 or i^o < ~1- 
The conditions for a gapped background solution are 

• Validity of the Gubser criterion < 3 ^ 

• Thermodynamical instability of the small black hole branch (^^ < 0): 5"^ > 1 

• The Sturm-Liouville problem is well-defined: |z^o| > 1 (stronger than It'ol > 1/2)- 

Figure 5 is a summary of the above discussion. The parameter space when the current- 
current correlator has a gap is region B. 
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4 Conductivity in the presence of momentum dissipation 


To add momentum dissipation and at the same time preserve the rotational symmetry we 
will consider the following bulk theory 


5 = y 


« 4- rw - ^ Ewy 


2=1 


, (4.1) 


where the leading IR behavior of V, Z, W, and Y are 

V{(f)) = Voe-^^, Z{(j)) = W{(j)) = Woe^^, Y{4>) = 


(4.2) 


This is a simple way to add momentum dissipation [45, 46] as the scalars V'i take the form 
V'i = kxi- Having the same slope for the axions guarantees a rotationally invariant solution. 
To calculate the conductivity, we perturb the system by 

5A^ = ax{r)e~"'^\ 6gtx = gtx{r)e~''‘^\ h'f/'i = y(r)e"*‘^^ (4.3) 


The independent linearized equations are 

< + ((logZ)' - ^ 4 + (5^ - Bax + § [g'tx - ^9tx ] = 0, 

C ik 

9tx - -;^9tx + A'^Zax H- DYx' = 0. 

O LU 


We dehne 


(jj \ B 


^/Im * 


(4.4) 

(4.5) 

(4.6) 

(4.7) 


From (A. 12), we observe that gr is a constant if VF = 0. If VF 0, g is the minus charge 
density in the UV, but not a constant in the bulk. In the fractionalized phase, g is a nonzero 
constant in the IR, while in the cohesive phase, g = 0 in the IR. After eliminating gtx{f)^ 
two coupled equations are derived 



ax - q- 


y/BD 


bx = 0 , 


-h'] _ 

CY\I B y ^ yCYy D j 


bx - k^q 


y/BD 

C2 


ax = 0 . 


(4.8) 

(4.9) 


The equations for ax and bx can be written in a more symmetric form. Dehne Zi := 
Z = and Z 2 := The relation between the metric (3.1) and (2.3) is 




















Also we define / = gfVh = \JDjB. The coupled equations for ax and hx 


are 






W 


Ux + qba 




Define Ux = \/~Ziax and bx = ^/Z^bx-, and then we obtain 


- 2~ , (igr 

- 7 :^ + bl(?)®a; — LO ttx -\ - 


de 

d?bx 


+ V2{^)bx — uj'^bx + 


h\/ Z1Z2 

j^2qgr2 

h\/ Z\ Z 2 


bx ) 


5 


where x is the same as in (3.5) and 

1 


^*(0 = 


i?2 


4 * 2 * 4 F2 2 F 




Pi = - 


{ fz^y 


Qi = ^ + 


Q 2 — 


k'ir^ 


yzr ^ gZi r^gZi' gZ2 

After inserting the hyperscaling-violating geometry to the equations, we obtain 


where 


d'^dx , c^„^ 

de [e ev 

d'^bx .(c2 C 3 \ 

de J 


~ — 2 ~ , di 7 

^ o,x T ^^bxi 

bx — ^ bx ~\~ ~^dxi 


(7 - (i )(7 + (5 + 2A) _ kX-2 
(7 - (5) ( 7 -I-35)-I-4 22 ; 

4(7 — 5)(7 — 5 — 7 ) K 7 -|-20 — 4 

(7 - 5) ( 7 -I-35)-f 4 z 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


(4.15) 


(4.16) 

(4.17) 

(4.18) 

(4.19) 


and Cl, C2, C3, c«,, di, and d 2 are constants.The coefficients ci and 02 are given by 
Cl = i/f — 1/4 and C 2 = t'l — 1/4, where 


(7 — 5) (37 -|- 55) -|- 12 3z — 9 
2 [(7 - <5) (7 + 35)-I-4] 22 ; 

(7 — 5) (7 — 55 — 4A) — 4 2 — z — 9 — k\ 

2[(7 - ^)(7 + 35)-I-4] 22 ; 


(4.20) 

(4.21) 


The terms involving a are from the momentum dissipation sector and the term involving 
(3 comes from the U(l) symmetry breaking. We impose two conditions: 

^^The coefficients C3, c™, di, and 6.2 are obtained from 


C3 

c» W 

di 

qgp 

d 2 

7 2 2 
k qgr 

^“-2 gZ 2 ’ 

PgZi ’ 


hy/Z^ZZ 


hx/ Zi Z2 


20 


































• The axions do not change the IR geometry to leading order. After we include the 
axions, the effective potential for the scalar is 

Hot = V{ct>) + ^g^^k^ ~ ~ ^-r5 as r ^ oo. (4.22) 

For the axions to not change the IR geometry, the first term must be dominant. 
Therefore we require 

— R5 > 2 + RX ^ (7 — ( 5)(7 + 5 + 2A) < 0. (4.23) 

• The system is in the fractionalized phase, in which the symmetry breaking term is 
subleading in the IR.^^ To implement this we demand that in (A. 12) the right hand 
side is subleading to the left hand side. This implies 

(7 — ( 5)(7 — 5 — 7 ) > 0. (4.24) 


If the IR is at ^ — 7 - 00 , i.e., the extremal limit is at T —)• 0, we have a > 2 and (3 > 2. 
If the IR is at ^ 0, i.e., the extremal limit is at T —)> 00 , we have a < 2 and /3 < 2. We 
will show that the terms involving a and (3 are all subleading to the solution without those 
terms. This will imply that the qualitative features of the charged spectrum is not changed 
by the presence of the momentum-dissipating sector nor by the U(l) symmetry breaking. 

We take first = 0 {W = 0) for simplicity. We can obtain a fourth-order equation for 
CLx after eliminating bx'- 


d^i 




+ 


2a d?d^ 


+ 


a{a — 1) — Cl — C2 ec3 \ d'^dx 2(2 — q:)ci dd. 


+ 


^3 ^2a-2 J ^^2 

Ci(c 2 - (2 - a)(3 - a)) 


+ 


+ 


^2 di 

e(ciC3 - did2) 


t20L 


dx = Q, (4.25) 


where e is a bookkeeping parameter that keeps track of the order of perturbation and which 
we will eventually set to 1. The asymptotic behavior of Oa, at ^ 0 is 


dx = (ao + eaie^( 2 -a) ^ g2^^^4(2-a) + ... ^ , (4,26) 


Since a < 2 , the above expansion is valid and reliable. The coefficients ai, 02 , • • • can 
be solved order by order. We observe that the coupled equations for dx and bx do not 
change the leading exponent of dx- Therefore, if dx is normalizable at ^ 0 without the 

momentum dissipation, it is still normalizable with the momentum dissipation. The hard 
gap remains in the conductivity. 

When the U(l) symmetry breaking term is included, the fourth-order equation for dx 
is 


2q 

+ + 


-t 


/ Oi(^Oi 1 ) Cl C2 €q-C3 

d^dx 1 

f'2(2 - a)ci 2ep{j3 - a)cw'\ 

ddx 

^2 ^2a-2 J 

' dH 

\ c ) 

• dc 


^ dC3 

ci(c2 - (2 - a)(3 - a)) ea(ciC3 - did2) e^c™(c2 - - o)(;0-t 1 - a)) tatpczCn 

j>-A ' >-0^, ' ' 


^/3+2 


^2a+/3- 


dx — 0 , 

(4.27) 


®In the cohesive phase, the symmetry breaking term is at the same order, which requires rj = j — 5. 
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where there are two subleading terms multiplied by eo = that in the end will be set to 
1. The solution is a double expansion 

Qx = ^ao + ^ ^ + ' ' ' ^ ) (4.28) 

where all terms involving a and /3 are subleading. Therefore, we confirm that in this case, 
the U(l) symmetry breaking term and the momentum dissipation effects do not charge the 
main characteristics of the spectrum: the spectrum remains discrete. The only difference is 
that the zero mode (responsible for the (5-function at zero frequency) has now disappeared 
and therefore now the spectrum is gapped. We conclude that the parameter space for the 
gapped geometry in figure 4, for the fractionalized phase has not changed. 

In the following we will discuss the solutions of ax and bx in the IR. We will impose 
the Gubser criterion so that the extremal limit is at the T —)• oo end of the small black-hole 
branch (^ —?■ 0). The general solutions of dx and hx are 

dx = ^ ^ (4.29) 

bx = . (4.30) 

For the system to be holographically reliable we must impose the condition that only one 
normalizable solution exists for each of dx and bx- This requires \ui\ > 1 and \u 2 \ > 1. We 
have the following observations: 

• < 0 is always satisfied, but i^i < — 1 gives nontrivial constraints for ( 7 ,( 5 ), as 
discussed before. 

• 1^2 can be either positive or negative, and |z^ 2 | > 1 gives nontrivial constrains for ( 7 , 6). 

In the UV, substituting in (2.16), the near-boundary expansion of the potentials in 
(4.15) is 

Vi = ^A(A - l)7(/>ar^"^ -F ^(A%^ -F 

+ ^(2 - A)(3 - A)j(j)br^-^ + ^((3 - A)V + 

V2=- ^A(A + l)A((.,r^-2 + ^A2A2(/,2r2(A-i) 

- ^(3 - A)(4 - A)X(j)br^-^ + ^(3 - 

The UV behavior of and bx is 

Ux = -I- -I-, bx = b^x'>r -\ -. (4.33) 

The Schrodinger potentials Vi and V 2 are divergent when 1/2 < A < 2, except for A = 1 
or 2, for which they asymptote to constants. The analysis is similar to the last paragraph 
of section 3.1. 


(4.31) 

(4.32) 


22 - 


5 On holographic Mott-like insulators and supersolids 


We will first consider the EMD bulk action with no U(l) symmetry breaking {W = 0 in 
(4.1)). To decide whether the system is a metal or insulator, we must calculate its DC 
conductivity. With momentum dissipation, a formula for the DC conductivity was derived 
in [42]. We will briefly review the procedure, and use it for the near-extremal geometries. 
Then we will apply this formulation to extremal geometries. We define for latter convenience 

//(r) = Z, + |lz,^Z+^. (5,1) 

The two eigenmodes of the equations for ax and hx are 


\ — ^ ( _ 3_^h 


Z 2 

^2 = ^ {qax + hx). 

tL 


(5.2) 


Following [42, 79], we can derive an equation 




OJ 


+ —H\i = 0. 


(5.3) 


From this we observe that the equations for ax and bx have a massless eigenmode. Therefore, 
we define 

n = fHX[ - ^fZ2 (A 2 , (5.4) 

and n is a radially conserved quantity at cj = 0. At the UV, we have H = 1,^® Ai = 

A 2 = 0, and n = a^^ in the leading order. The conductivity is 


a{u;) 


ax'* 

■ ( 0 ) 
iLoax 


n 


iuXi 


r—>0 


At w = 0, n can be evaluated at the horizon. 


(5.5) 


<7dC 


n 


itaXi 


r=rh 


(5.6) 


Refs. [79, 80] derived the DC conductivity at finite temperature as 


c^DC = Zh + 


q 


k^ChYh 


= + 


q 




(5.7) 


where (l)h = and Ch = C^r^) are the values of the scalar (dilaton) and the scale factor 

on the black-hole horizon. At finite temperature, ctdc is finite. The formula in (5.7) however 
cannot be directly applied to extremal geometries. We will discuss the finite temperature, 
near-extremal geometries first, and then the extremal geometries. 

^®As we set Zq — 1. 
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5.1 DC conductivity for near-extremal black holes 


We now return to the analysis of the nature of the conductivity in the hyperscaling violating 
Lifshitz geometries perturbed by U(l) breaking as well as by the momentum dissipating 
sector. 

We apply the formula (5.7) for the DC conductivity to the near-extremal geometries, 
and then take the extremal limit. The near-extremal solution from [17] is^® 

ds^ = —V{f)f dt^ -t- 


-.0 ,0 7 2 \ 

+ r {dx +dy), 


= e‘^°f ^ 


( 7 -^) 


—wAV{r) 

At = 2\ -fh) V{f) = r{r - rt), 

V wu 

u = — ■yd + 2, V = 5^ — 'yd — 2. 


wu = 37 ^ — 2yd — <5^ -|- 4, 

Note that this r is different from the r in (2.3). Thus we obtain 

.,2 _ 2(-j-S)(j-S-2\) 2 

- 3^2-2-15-52+4. _ ~Ol i ^ 

p'h - ' h ^ J^2‘ h 

In terms of z and 9, 


— 47(7-6) 2 _ 2 (-/- 5 )(-f- 5 - 2 \) 

n-T^x^ r-u r 37^-276-6^+4 Wf 373-276-63+4 _ ~ai I H -012 

ODC ' h — I h > ~ ^ ' 1. 


e-A 


Oti — 


0-2 = 


9 -2 +k\ 


(5.10) 


(5.11) 


2 + Z - 9 ' " 2 + Z -9 

The extremal limit is r/j —)• 0 keeping the other coordinates hxed. The Gubser criterion and 
the condition that the axions are irrelevant in the IR imply that the second term is always 
dominant: 

2 ( 7-<^)(7 + '^ + 2 A ) 

372 -2'yd -d^ + A ^ ' 

Therefore, the DC conductivity is always dissipation-dominated in the extremal limit. If 
(Tdc —)• 0 as f/j —>■ 0, the system is an insulator; if ctdc 00 , the system is a (perfect) 
conductor. The condition that the system is an insulator is 02 < 0, i.e.. 


(7 — h) (7 — (5 — 2A) < 0. 


(5.13) 


For different values of the parameter A that controls the IR strength of momentum dissi¬ 
pation, the parameter spaces for insulators and conductors are plotted in hgure 6. In the 
green region, the system is a conductor, and in the red region, the system is an insulator. 
Finally in the light blue region the axions are relevant in the IR, a case that was analyzed 
further in [79, 80], but we do not analyze it further here. 

The relation between the temperature T and the horizon radius r/j is [17] 


T 


1 

Att 


V—wAe 


( 7 - 7 ^ 


7^+276-363+4 

-3+3-276-63+4 


(5.14) 


^°The AdS 4 near-extremal solution in coordinates (2.3) is (we need to add appropriate coefficients before 
dt^ and dr^, as in eq. (2.4)) 


ds^ = 


1 


/ 4(z-l) 25 

l-fr 2 -b dt Jr r 2 -s 


dr^ 

T 



/ = !- 



2(2 + z-6) 
2-6 


(5.8) 


The power in the blackness function /(r) is always positive as required by the Gubser criterion. The 
extremal solution can be obtained by sending r/i —>■ 00 . 
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Thus, the DC conductivity for the near-extremal black hole is 


_ 47(7-(5) _ 2(7-(5)(7-5-2A) o 

<7DC ~ T 7^+275-35^+4 _)_ 72+27^-352+4 = J'Pl _|_ ^ 


In terms of z and 6, 


/3i = 


A:2' 

6»-4 C-2 




/32 = 


0 - 2 -h kA 


(5.15) 


(5.16) 


where C is the conduction exponent given in (2.11), [19, 71]. As mentioned above, the 
second term is always dominant no matter whether the extremal limit is T —)• 0 or T —?■ oo. 

We are especially interested in Mott-like insulators. In figure 6 , the parameter space 
in which the conductivity can have a hard gap and a discrete spectrum is above and below 
the two curves (7 — h )(7 -|- 3h) -|- 4 = 0 (which gives = 0). This is in the region B in 
figure 4. From figure 6 , we can see that this region can be either a perfect conductor or an 
insulator, depending on the value of A. In sum, in the red region of the parameter space for 
the gapped geometry, the system is a Mott insulator, if the scaling dimension of the scalar 
operator in the UV satishes 1/2 < A < 2. 

If the U(l) symmetry is spontaneously broken, there will be a h-function at a; = 0 
even in the presence of momentum dissipation (the superfluid mode) provided that the 
breaking is spontaneous and we will have therefore a superfluid (or superconductor if the 
U(l) is weakly gauged) with a discrete spectrum. This will provide a system that it is 
in the broken phase there will be interpreted as the modulus of a charged scalar), 
and will therefore have a (5-function at zero frequency in the conductivity but will have 
discrete spectrum otherwise! This characteristic ressembles what is expected from the 
conductivity in supersolids (although there are also differences that were expanded upon in 
the introductory sections). 


5.2 DC conductivity for gapped, extremal geometries 

For extremal geometries, the formula for the DC conductivity (5.7) cannot be directly 
applied. When the extremal limit is at T —>■ 0, the terms in the perturbation equations 
can not be neglected to calculate the cu —)• 0 behavior of the conductivity; we need to use 
the asymptotic match. On the contrary, when the extremal limit is at T —)> 00 , we can 
neglect the terms. Moreover, for gapped geometries, the IR solution is independent of 
LO, and there is no in-falling wave boundary condition. From (5.5), we can evaluate the 
radially conserved quantity 11 at the IR, and the a; —)• 0 limit of the conductivity is 


croc 



iojX 


( 0 ) ’ 
1 


(5.17) 


where A® = is the source of the perturbation, and 11 can be regarded as the strength 
(coefficient) of the zero-frequency (5-function. 

The condition that the presence of axion-generated dissipation do not change the orig¬ 
inal IR solution is equivalent to 


Z1/Z2 —)• 0 


as 


r —)• 00 . 


(5.18) 
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Figure 6. Parameter space for the conducting and insulating phases. The parameter space 
constrained by the Gubser criterion is enclosed by dashed lines; moreover, we have excluded the 
holographically unreliable region in which there are two normalizable solutions in the IR. In the 
gray region, the axions term is relevant in the IR. In the green region, the system is a conductor; in 
the red region, the system is an insulator. The parameter space for the gapped geometry is above 
and below the two curves (7 — i 5)(7 + 3(5) + 4 = 0. 


Another important observation is that we always have^^ 


Z\ (lx 
Z2 bx 


r\j 



^Wl\-W2\ ^ 0 , 


(5.19) 


^We do not always have a^/bx ^ 0 as ^ 


0 . 
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no matter 1/2 > 0 or 2/2 < 0. Consequently, the leading behavior of 11 in the IR is 


Il = -^fZ2b', + --- (5.20) 

If n —7- 0 as r —)• 00 , the conductivity does not have a d-function at w = 0. If 11 is constant 
as r —)• 00 , the conductivity has a d-function at a; = 0 and its coefficient is proportional to 
n(a; = 0). That there is only one normalizable solution for bx requires 11 ^ 2 ! > 1- There are 
two cases to consider as follows: 

• Case 1: 1/2 > 1, for which we have (7 — S){'y — 55 — 4A) — 4 < 0. The solution for bx 

is bx ~ The radially conserved quantity 11 evaluated at the IR is 

(7 —5) (7—55—4 A) —4 

n ~ fZ 2 b'x ~ r (TTp _i. 0 . (5.21) 

• Case 2: 1/2 < —1, for which we have (7 — 5)('y — 55 — 4A) — 4 > 0. The solution for bx 

is bx ~ . The radially conserved quantity evaluated in the IR is 

n ~ fZ 2 b'x —)• constant. (5.22) 


We expect that the insulating phase corresponds to the case 1, and the conducting 
phase corresponds to the case 2. We will show that this is true after we excluded the 
holographically unacceptable region \i/ 2 \ < 1 and have imposed Gubser’s criterion. 

Gapped, insulating phase. The condition that the extremal geometries are at T —)■ 00 
limit of the small black-hole branch is (7 —5) ( 7 -|-35)-|-4 < 0. The condition that the system 
is an insulator is (j—5)(j—5—2A) < 0. These two inequalities give ( 7 — 5 )( 37 -|- 5 —4A)-|-4 < 0. 
Consequently, 


2^2 + 1 


(7 — 5) (37 -I- 5 — 4A) -I- 4 
(7 - 5)(7 -I- 35) -I- 4 


(5.23) 


If we further impose that there is only one normalizable solution, which is \i' 2 \ > Ij we 
obtain 1/2 > 1- According to the discussion above, this corresponds to the case 1. The 
conductivity does not have a 5-function at w = 0 in the gapped, insulating phase. Note 
that this is not always true without the constraint |jz 2 | > 1 . 

Gapped, condueting phase. The condition that the extremal geometries are at T —)■ 00 is 
( 7 — 5 )( 7 -|- 35)-|-4 < 0. The condition that the system is a conductor is ( 7 — 5 )( 7 —5—2A) > 0. 
These two inequalities give (7 — 5) (7 -|- 115 -|- 4A) -|- 12 < 0. Consequently, 


(7-5)(7 + 115 + 4A) + 12 
(7 - 5) (7 -I- 35) -I- 4 


(5.24) 


If we further impose that there is only one normalizable solution, which is \u 2 \ > 1, we 
obtain U 2 < —1. According to the discussion above, this corresponds to the case 2. The 
conductivity has a 5-function at a; = 0 in the gapped conducting phase. 
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6 Constructing the finite temperature geometry 

To determine which is the ground state of the system (2.1), we need to solve for the hnite 
temperature geometry, and lower the temperature to obtain the near-extremal geometry. 
Note that the RN black hole is not a solution of the system (2.1) when Z'((/>)|^=o 7^ 0- 
If we want to choose A = 1 for general values of 6 without log r terms in the boundary, 
we can use a three-exponential potential. The simplest potential is 

R((/>) =-6-psinh^ y. (6.1) 

The reason why we choose this potential is explained in appendix B. 

After eliminating A'/- by '/hZ{cj))A[ = —p, where p is the charge density, we solve the 
three coupled equations for cj), g, and h. We work in the canonical ensemble. At hnite 
temperature, we can set the horizon be at r/i = 1 by a rescaling of coordinates. The 
asymptotic behavior of the scalar held (p near the AdS boundary is 


(j) = (j)ar + Pbr'^ H-. (6.2) 

We consider the UV scaling dimension A = 1, so is the expectation value, and pb is the 
source. The way to impose the boundary condition in the UV is not unique, for example, 

• (pb = (po- This is the standard Dirichlet case. Here po has the dimension of energy. 
We have a two-parameter family of the solutions (p, po). 

• Pb/Pa = 1^- This corresponds to a double-trace perturbation by the scalar operator 
dual to p. In this case, k has the dimension of energy. We have a two-parameter 
family of the solutions (p, k). 

• pb/Pl, = T- This corresponds to a triple-trace perturbation by the scalar operator 
dual to p. In this case, the only quantity that has a dimension is the charge density. 
We have a two-parameter family of the solutions (p, r). 

On the scalar we will use boundary conditions of the third type.^^ We further need four 
boundary conditions as in the following table (“X" means no boundary condition). 



horizon (r = 1) 

AdS boundary (r = 0) 

p 

regularity 

p" - 2 tP'^ = 0 (A = 1) 

At 

At = 0 

A't = -P 

9 

9 = 0 

X 

h 

X 

h = l 


^^The bulk action that has the gapped solution, has also an AdS2 IR critical solution. It is RG-unstable. 
The bulk action that has the gapless solution has also an AdS2 IR critical solution which is RG-stable. 
However we do not reach it for the values of the sources we explored in our numerics. We should also 
note that the parameters of the action with the gapped solution were chosen slightly outside the region 
determined in the previous section for reasons of numerical reliability. 
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Figure 7. Numerical results for physical properties of a gapless geometry. The parameters are 
chosen by 7 = 0.5 and 6 = —0.7. In this case, the extremal limit is at T —0. 


The entropy is S' = = A-kV, as we have set rh = 1, where V is the volume (area) 

in the spatial dimensions. The dimensionless temperature and entropy density are 


f = — = 



(6.3) 


The extremal limit is T —)• 0 or T —)• 00 . Although the RN black hole has finite entropy at 
zero temperature, the other cases studied here have S —)• 0 in the extremal limit, which is 
obtained by sending p —>■ 00 . When the geometry is near-extremal, we obtain 

S(xf^. (6.4) 


The free energy in the canonical ensemble is 



^ T 

(6.5) 

where e is the energy density from 

9 

h 

= 1 - -r^ H-. 

r^O 2 

( 6 . 6 ) 


We compare a gapless case and a gapped case. We choose the parameters ( 7 , 5) ac¬ 
cording to figure 4, and solve the equations of motion numerically. 

Figure 7 is for a gapless case with ( 7 , 6) = (0.5, —0.7), and figure 8 is for a gapped case 
with ( 7 , 5) = (0.4, —1.2). In figures 7 and 8 , the upper-left plot shows the relation between 

















Figure 8. Numerical results for physical properties of a gapped geometry. The parameters are 
chosen by 7 = 0.4 and 6 = —1.2. In this case, the extremal limit is at T — 00 . 




Figure 9. AC Conductivity calculated at two different temperatures: = 0.538 for the left 

plot, and = 0.858 for the right plot. The background is a gapped geometry with 7 = 0.4 

and S = — 1 . 2 . As we go to the extremal limit by increasing the temperature of the unstable branch, 
the spectrum will have a hard gap. 


the temperature and = 4>{rh)- The upper-right plot shows the relation between the 
entropy and the temperature. The lower-left plot shows the scaling exponent in S ~ as 
a function of (ph, and the lower-right plot shows the relation between the free energy and 
the temperature. 

In the gapped case, there is a minimal temperature T^m for the black holes and there is 
a small black hole and large black hole branch. This is clearly visible in the left-upper plot 
of hgure 8 . The small black hole branch is on the right side of this hgure. The extremal limit 
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can be approached from this branch by going to the right. In this limit the temperature 
of the black holes T —)• oo like the usual Schwarzschild black holes. Only the left branch is 
stable. 

The extremal limit is the ground state of the system. As we heat up the system, there 
will be a first-order phase transition at Tc > Tmin and the system will jump from the 
low temperature ground state to the black hole in the stable branch, [56, 57]. The AC 
conductivity calculated from the same gapped geometry at two different temperatures is 
plotted in figure 9. It shows that as we increase the temperature of the unstable branch in 
order to approach the extremal limit, a hard gap will emerge in the extremal limit. 


7 The spectrum of tensor and scalar fluctuations 


In previous sections, we studied systems where the current-current correlator has a gap. 
Here we study the spin-2 and spin-0 excitations associated to the stress tensor. Such 
fluctuations determine holographically the spin-2 (also known as the shear channel) and 
the spin-0 (bulk channel) parts of the stress-tensor two-point function. 

Generally, we need to turn on the following perturbations and derive the linearized 
equations of motion, to calculate various transport coefficients: 

(1-|-e *‘^*/i^j/(r)), —)• A^(l-|-e *^*a^(r)), (/>—)> 0 (1-|-e 

(7.1) 

The fluctuations can be classified by their spins under the spatial rotations (here i = x, y) 
as follows. 

• spin-0: hu, h^r, ha, htr, at, ar, 6(f>. In a gauge where 5(j) is unperturbed [81, 82], a 
decoupled equation for ha can be obtained, responsible for the bulk viscosity. 

• spin-1: hti, hri, a*. A decoupled equation for at can be obtained, responsible for the 
conductivity. 

• spin-2: hij. The equation for hij is responsible for the shear viscosity. 


7.1 Spin-2 fluctuation - Shear viscosity 

The spin-2 fluctuation hxy satisfies the Laplace’s equation. The Laplace’s equation V^<I> = 0 
with the the metric (3.1) is 


+Y^Ct<;2$ = 0. 


After the change of variables by 

di _ [B 

dr y D’ 

we can obtain a Schrodinger equation 

d^^ 


= Vc^, 


de 


+ H(^)$ = 


(7.2) 


(7.3) 


(7.4) 
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The potential is 


DB'C D'C DC'"^ DC” 
4B^C ^ ABC ~ ABC^ ^ 2BC' 


(7.5) 


where the prime is with respect to r. If V goes to to infinity in both IR and UV, then the 
Schrodinger eqnation can only have bound states that correspond to a discrete spectrum 
for the spectral density of the correlator. To analyze the behavior of V in the IR, we use 
the hyperscaling violating geometry (2.4) as the background. 

IR charged solution. We change the variable r to ^ by (3.9), and the Schrodinger 
potential is 


V{x) = 



(7.6) 


where 

_ (7-(5)(37 + (I) + 4 _ 2 + z - 9 

2[(7 - h )(7 + 3h) + 4] 2z 


(7.7) 


Assuming the Gubser criterion is always satisfied, we have vl > 0 when the extremal limit 
at T —)• 0 (^ —)• oo), and < 0 when the extremal limit is at T —)• oo (.^ —)• 0). In the later 
case, we should further impose that |i/| > 1 so that there is only one normalizable solution 
in the IR. By comparing (3.11) with (7.7), we can see that \v\ — \vl\ > 0, which implies 
that when the Laplacian is gapped, the current-current correlator is also gapped and vice 
versa in the holographically reliable region. Figure 10 is a summary of the parameter space 
( 7 , 6) for the spectrum of the Laplacian. 

IR neutral solution. We change the variable r to ^ by (3.9), and the Schrodinger 
potential is 


V{x) 


^2 ’ 


(7.8) 


where 


h2-3 _ 3-9 

^^0 “ 2 { S 2 _ 1) - 2 


(7.9) 


We observe that no in terms of 9 can be obtained from (7.7) with z = 1, although the 
expression in terms of 6 is very different. Assuming the Gubser criterion is always satished, 
we have n^o > 0 when the extremal limit at T —>■ 0 (^ —)> 00 ), and nio < 0 when the 
extremal limit is at T —)> 00 (.^ —)• 0). In the later case, we should further impose that 
Wlo\ > 1 so that there is only one normalizable solution in the IR. Gonsidering all above 
constraints, the parameter space for the gapped Laplacian is 1 < (5^ < 5/3. 


7.2 Spin-0 fluctuation - Bulk viscosity 

To calculate the bulk viscosity, we can turn on the perturbations hu, hrr, hxx, and at only, 
according to the gauge choice of [81, 82]. After eliminating hu, hrr, and at, a decoupled 
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Figure 10. The regions A, B, and C are the parameter space constrained by the Gubser criterion. 
We assume that the conditions in the UV for a gapped geometry has been satisfied as explained in 
the end of section 3.1. In region A (yellow), the extremal limit is at T —>■ 0, and the Laplacian is 
gapless. In region B (red), the extremal limit is at T —>■ oo, and the Laplacian is gapped. In region 
C (green), the extremal limit is at T —)■ oo, and the Laplacian is gapless. Region D (enclosed by 
blue boundaries) is holographically unreliable. When both correlators are holographically reliable 
a gapped current-current correlator implies a gapped shear tensor correlator. 


equation for h^x is derived^^ 
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Kx = ^- (7.11) 


di 

dr 


y - \ 1 dxX - ^ hx 


where the prime is with respect to r, we can obtain a Schrodinger equation 

d'^hx 


d^ 


T ^{^)hxx — ^ hx 


(7.12) 


(7.13) 


If V goes to to infinity in both IR and UV, then the Schrodinger equation can only have 
bound states that correspond to a discrete spectrum for the spectral density of the corre- 

^^This appearance of this equation is not unique by applying the equations of motion. It can be written 
in a simple form in the coordinates (2.3) with / = g/Vh: 


1 

w 





r^hjhZyA'y 

gh' 


jg/h)'{r<j>'y \ 
g/h rtj)' j 


hxx — 0. 


(7.10) 
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lator. To analyze the behavior of V in the IR, we use the hyperscaling violating geometry 
(2.4) as the background. 

m charged solution. We change the variable r to ^ by (7.12), and the leading order of 
the Schrodinger potential is 

— 1/4 

^(0 = ( 7 - 14 ) 

where 

_ V[(7-'5)(37 + <5) + 4][(7-<5)(37 + 175) + 36] _ ^/(z - 0 + 2)(9z - 9 - 6) 

""" 2[(7-,5)(7 + 35)+4] 2zsgn(2-0) ’ 

Assuming the Gubser criterion is always satisfied, the expression inside the square root is 
always positive, and we have Ug > 0 when the extremal limit at T —>■ 0 (^ —)> oo), and 
Ug < 0 when the extremal limit is at T —>■ oo (^ —)• 0). In the later case, we should further 
impose that > 1 so that there is only one normalizable solution in the IR. The three 
correlators responsible for the conductivity, shear viscosity, and bulk viscosity are gapped 
at the same time in their common holographically reliable region. 

IR neutral solution. The Schrodinger potential is the same as (B.7). Therefore, the 
discussion about the spectrum for the scalar fluctuation is the same as the Laplacian. 

8 Conclusion and outlook 

We have obtained two new classes of holographic matter, constructed from geometries that 
have correlators with discrete spectrum of excitations. 

The first is a class of holographic ground-state saddle points of the bulk action (1.3) 
with W{4>) = 0 which implies the U(l) symmetry is intact. Such saddle points have a source 
of momentum dissipation due to the non-trivial axion solutions that break the translational 
symmetry. The characteristic exponent A in (4.2) that controls the IR properties of the 
coupling of the axions to the rest of the degrees of freedom can be selected so that the 
momentum dissipation effects are irrelevant in the IR. 

In the absence of momentum dissipation, the saddle-point solutions of interest, found 
already in [17] have a discrete spectrum for the charge charge correlators, implying the 
associated spectral density is a discrete sum of (5-functions. Translation invariance implies 
that one of the (5-functions is at zero frequency. 

A careful analysis, done in this paper, indicates that in the presence of IR irrelevant 
momentum dissipation, the correlator still has a discrete spectrum, with the only difference 
that the zero frequency 5-function has now disappeared. Therefore, these are bona-fide 
insulators with a hard gap and a discrete spectrum of charged degrees of freedom. They 
resemble in several respects Mott insulators but we suspect this is a novel class of insulators. 

We have found explicit numerical solutions that connect well-defined asymptotically 
AdS 4 UV fixed points with a non-trivial charge density to the insulating fixed points dis¬ 
cussed above. There are however caveats in our discussion. Although explicit flows exist, it 
is not yet clear whether they are the dominant (in the sense of free energy) flows between 
the UV and IR fixed points in top-down string theory effective actions, although the M- 
theory example of [55] is a good candidate for a stable flow. We believe that the mechanism 
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is a wider interest and is the one expected to be at play for example in real QCD at finite 
but 0(1) baryon density, [58, 83]. 

Other characteristics of such saddle points are as follows: 

• All non-trivial quantum critical saddle points with z ^ 1 have a mild IR geometrical 
singularity. It is a naked singularity in the presence of hyperscaling violation, while 
it is a much subtler singularity not captured by scalar curvature invariant in the 
case of Lifshitz scaling symmetry. Despite the fact that holography demands regular 
saddle points, singularities can be tolerated because they can arise from the neglect 
of some degrees of freedom (coarse graining). They should be however resolvable, as 
adding back the missing degrees of freedom, one should recover a regular background. 
A criterion for resolvability was advanced by Gubser, [76], and we always impose 
it on our singular saddle points. There is however a further issue in the presence 
of a resolvable singularity: that the calculation of correlators should not depend on 
the resolution of the IR singularity. This issue crops up as the presence (or not) 
of two normalizable modes, in the Sturm-Liouville problem that is relevant for the 
calculation of a two-point function. In such a case the solution is not unique, and 
a further boundary condition at the singularity is necessary in order to define the 
correlator, [17]. This is the case where the holographic calculation is unreliable, and 
a resolution of the singularity is necessary in order to provide a concrete answer for 
the relevant correlator. 

We find that in the saddle-points of interest studied in this paper, the (5-function at 
zero frequency, is absent only for those bulk theories for which the charge correlator 
is holographically well defined (does not suffer from singularity resolution ambigu¬ 
ities). This is interesting and may provide important clues for other cases where 
zero-frequency (5-functions in transport appeared when it was not expected, [9]. 

• In a related theme, in appendix C we have computed analytically the conductivity of 
a class of black holes solutions of gauged supergravity (known as STU black holes) 
dual to deformations of the ABJM theory. In particular, in the extremal two-charge 
case (where the charge density is zero), we obtain a gap and a zero-frequency 5- 
function although none was expected due to symmetry principles. The reason is that 
there is no charge density and without charge density there is no argument for a 
zero-frequency (5-function. Moreover, the equation for the conductivity below the gap 
is in the holographically non-well-defined class as it has two normalizable solutions. 
However here a unique solution can be chosen using analyticity in the frequency oj. 

• The zero-temperature phase of the gapped solutions persists up to a finite critical 
temperature Tc- Indeed the phase diagram of candidate black-hole ground states is 
similar to that relevant for holographic QCD, [56-58]. There is a minimum temper¬ 
ature, Train < Tc below which there are no black-hole solutions as can be seen in 
the upper left part of figure 8. Above Tc the system jumps to a plasma-like phase 
via a first order transition. Calculating the conductivity of this phase is an interest- 
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ing problem. We expect that this will correspond to a conductor, where momentum 
dissipation effects may be strong. 

• In the low temperature insulating phase the change of conductivity with temperature 
is an effect suppressed in the large-limit. It can be calculated however by a one-loop 
gravity calculation. The relevant diagrams will be the one-loop correction to the bulk 
propagator of two gauge bosons. This calculation is in principle feasible (a related 
one-loop correction to holographic thermodynamics was performed recently in [84]). 

• Energy transport is controlled by the stress tensor two-point function. In section 
7 we analyzed the equations for shear and bulk channels. We found that for all 
holographically reliable geometries, when the current-current correlator is gapped, 
then the shear and bulk channels are also gapped. Therefore there is no low-energy 
transport in these insulators. On the other hand, in the supersolid class this fact is 
consistent with our interpretation. 

• Our analysis has used effective holographic theories following the philosophy devel¬ 
oped in [17]. Some of the solutions we find and analyze however appear in string- 
derived supergravities, (see [55] and appendix C), namely the AdS 4 gauged supergrav¬ 
ity obtained from compactification of M-theory on 5^ that is dual to the 3d ABJM 
theory, [85] . It is therefore expected that they may be ubiquitous in three dimensional 
QFTs that are deformations of the ABJM theory at finite density. 

We do expect that in codimension-one subspaces of the theory space we are studying 
there will be IR fixed points that are non-scaling. This is the case we expect to give rise to 
Arhenius-type insulators and looking for them is an interesting endeavour. 

The second class of holographic saddle-points have finite density but in theories where 
the U(l) symmetry is broken. The action in (1.3) with W{(p) 7 ^ 0 is describing this class of 
saddle points and the universality classes of quantum critical IR fixed points were classified 
in [19]. This class contains both superffuid IR fixed points when the U(l) symmetry is 
spontaneously broken as well as fixed points with the U(l) symmetry explicitly broken. 

There are several possible IR fixed points but in this paper we have studied the frac¬ 
tionalized phases for which the bulk mass term W{4>) flows to zero in the IR. There is still 
a source of breaking of translational invariance in the form of non-trivial axion profiles. 

We have found fractionalized fixed points with a current-current correlator that has 
a discrete spectrum. Unlike the previous case studied above and despite the presence of 
momentum dissipation, the presence of symmetry breaking guarantees that there is a non¬ 
trivial (5-function at zero frequency due to the superffuid mode. 

Such systems are therefore superffuids (or superconductors when the U(l) is weakly 
gauged) in the presence of translational symmetry breaking with a discrete spectrum. This 
is not unlike to what is expected in supersolids with the sole exception that the translational 
invariance breaking is spontaneous. Like the case of insulators mentioned above, the systems 
here are also separated by a first order phase transition from a liquid-like deconfined plasma 
phase. Moreover the temperature dependence of the low-temperature supersolid-like phase 
is again a subleading effect in 1/N. 
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There are several problems that remain open in this direction. We highlight that the 
fate of other phases, in particular the cohesive ones requires further study. 

We also expect that the nature of the breaking of translational symmetry in the super¬ 
solid like ground states can be improved to be as required for standard supersolids, [59]. 
This can be done by “injecting" skyrmion-like charged solitons in the geometry (these are 
baryons in holographic models of QCD that are represented by flavor gauge theory instan- 
tons, [86, 87]). Such charge is expected to form a crystal at sufficiently high density, which 
is not expected to be high enough so that there is serious backreaction on the original 
geometry. Implementing this construction is an interesting problem. 

In this respect, looking back at standard (massless) QCD at zero density, we can 
interpret it as a zero density superfluid, where the broken symmetry is the SU{Nf)A flavor 
symmetry and the superfluid is composed of the massless pions. Turning on a baryon 
density, the system is driven to a quantum critical AdS 2 phase, [58]. It is plausible that for 
higher density there may be a crystalline phase that would be a supersolid. 
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A Equations of motion 


In this appendix we make explicit the equations of motion from the action we are using in 
(1.3) 


R,u - \g,.R = ^ , 


{Z{ct>)F^,) = W{ct>)A,, 


72 . W'icP)^,_Y'{<P) 


V^^iY{cj))V^A) = 0, 


-F^ - 


= 0 , 


2=1 


where 


= Z{cj>) (^F^pFP, - + W{cP) > 

Tt = d^(l)du(t) - \g^Jiu{^(|)f - gf,uV{(t>), 


2=1 


1 




(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 


We take tpi = kiXi {i = 1,2), which solves (A.4). The equations of motion with the ansatz 
ds^ = —D{r)dt^ + B{r)dr‘^ + C{r){dx‘^ + dy'^) (A.8) 

are as follows. The Einstein equations are 


2 ( 7 " a fB' a 
~Cr ~ ~C \ B ^ ^ 2D 


ZA'? BWAf PBY 1 ,,2 

+ ^ + + / =0’ 


2T»" 2(7" 


D C 
The Maxwell equation is 

The equation for the scalar is 


2(7" 
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(B' 
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(7 
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+ C 

D' 

(B' 
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D' 

D 

V B 

(7 

^ D 


BWA^ -2 

—^ = 0 , 


B'C 2ZA'? 2k‘^BY 


BC 


D 


C 


= 0 , 


f CZ 


B 


A'A =CWx -At, 


D 


+ 


C D' B'\ A',^Z'{4>) , AjBW'{ct>) k‘^BV'{cl>) 

C + + -~ 20 


Without the momentum dissipation (fc = 0), there is a conserved charge 

With the momentum dissipation, Q satisfies the following equation: 

Q! + k'^y/lmY = 0. 


(A.9) 
(A. 10) 
(ATI) 

(A.12) 

= 0. (A.13) 

(A. 14) 


(A.15) 
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B The potential for the scalar 


In this appendix we analyze optimal choices for the scalar potential. 

We assume that the potential for the scalar 1^(<^) only contains exponential terms, as 
is the case in supergravity systems. We need only the leading exponential term in the IR 
to discuss the hyperscaling violating geometry. However, we need at least two exponential 
terms to construct an asymptotically AdS geometry. A potential with only two exponentials 
can be parameterized as 


1^-exp — 


6 

L2[6J2 + A(3-A)] 


A(3 - A)e"'^('^"'^°) + 


(B.l) 


We have traded the four real parameters of a two-exponential potential in terms of the 
position of the UV fixed point, 4>o, the UV AdS^ scale L, the UV dimension A of the scalar 
operator dual to ip, and the extra exponent S that controls the IR hyperscaling violating 
solution when the factor dominates in the IR. 

The potential (B.l) (at zero charge density), has a single UV fixed point at cp = cpQ. 
In the IR, at finite charge density, it has two hyperscaling-violating fixed points. One 
corresponds to —5<p —>■ oo with an IR asymptotic geometry given in [17] with exponents 
( 7 ,( 5 ). The other corresponds to a hyperscaling violating geometry when 6(p ^ oo with 
exponents ( 7 , — 

This parametrization emerges from imposing the following constraints on the potential 
• U( 0 ) = —, Por A(i 54 with L = 1 , U( 0 ) = — 6 . 


• U'(0) = 0. 


• l/"( 0 ) = m2. For A = 1 or 2, U"(0) = -2. 


By applying the above three constraints, the potential only has one parameter that is 5. If 
we use the following potential 

y{4>) = (B. 2 ) 

which has V = —6 — cp"^ + ■ ■ ■ in the UV limit, the asymptotic behavior of (p will contain 
In 2 ; terms in general. The near horizon behavior of ip, g, and h is 

(p = hr + hh + hh’ H-h logr(/i2r2 -h (B.3) 


If we choose A = 2, there are no In 2 ; terms. If A = 1, the consistency of the equations 
without In 2 ; requires 5 = Tl/v/S. However, these values of 5 cannot give the gapped 
geometry. 

Consider a three-exponential potential 

V{(P) = cie"2<> + + cse^s'^. (B.4) 
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If we set cg = 0 for simplicity, we have one more free parameters in addition to <5. We 
parametrize the potential as 




2(6 + ri)e (6 + 

2 + (6 + u)5‘^ 2 + (6 + u)5‘^ 


(B.5) 


With this potential, if we plugin the near boundary expansion (B.3) without the log terms 
to the equations of motion, the leading term is 


2/i (~^ + (6 + u)6‘^) 2 

(6 + u)6 


(B. 6 ) 


If the near boundary expansion (B.3) is consistent without the log terms, the coefficient in 
the above term must be zero, which gives 


2(1 - 3 ( 12 ) 


(B.7) 


Substituting the above u to (B.5), we obtain the potential (6.1) used in our numerical 
calculations. 

Note that there may also be an AdS 2 x fixed point at constant cj) = (j)* [19], as 

shown in figure 1. The condition for the existence of an IR AdS 2 extremum at finite charge 
density is 



V 


(B.8) 


The critical exponents for the perturbation around the IR AdS 2 is 


P± = v^l - 4<5i^ , 


(B.9) 


where 


<li 


\v z vy 



(B.IO) 


If a RG-stable AdS 2 fixed point exists, the solution of (/>* is real, and there is an irrelevant 
mode associated with this (f)^. For the potential (6.1) and Z = the parameter space of 
( 7 , 6) under the above condition is potted in figure 11 . 


C Analytic solutions from supergravity 

Four special cases of the maximal gauged supergravity in AdS^ are given in Table 1, re¬ 
garding how they fit into the general classification. These four examples provide both the 
analytic understanding of some EMD systems and the numerical check of the program. The 
metric is 

^ + dy?) + ^2A{f)h{r) ' 

The relation between f and the r in (2.3) is = l/r2. The functions A{f) and h{f) are 

given in Table 1. These solutions are obtained from a more general black hole with four 
U(l) charges found in [ 88 ]. 


40 












2 

1 

5 0 

-1 

-2 

- 2-10 1 2 
K 

Figure 11. The parameter space of ( 7 , S) under the two conditions: (1) The solution of (j)* exists; 
(2) (3- < 0. In the light blue region, there is only one solution of (/)*, and in the blue region, there 
are two solutions of 0*. The parameter space depends on the concrete form of V[(j)) and and 

this plot is for the potential (6.1) and Z = 



The temperature and entropy are given by 

T = s = (C.2) 

dvr 

The relation between the free energy and the temperature is plotted in figure 12. There are 
some exceptional features for these systems. For example, if we set the horizon size rjj = 0, 
then the gauge field will vanish for the 1 -charge and 2-charge black holes. Consequently, the 
extremal geometries have to be at zero density. The geometries are significantly simplified, 
and it allows us to obtain exact, analytic solutions for (j{uj) for arbitrary a;. The extremal 
geometries for the 1 -charge and 2 -charge black holes can only be at zero density, and thus 
they are not “topologically" connected to the finite temperature black holes. The situation 
is similar to the asymptotically-AdSs case that was discussed in [89]. 

2-charge black hole in AdS^. To obtain the conductivity, we perturb the system ( 2 . 1 ) 
around the solution to (2.5) by 5Ax = ax(f)e“*‘^*, and obtain the equation for ax- 


ax + 


i , UJ 

f -b f2(f -b <5)2 


Qx — 0 , 


(C.3) 


where the prime stands for derivative with respect to f. By the usual redefinitions this 
equation can be mapped to a Schrodinger equation with a constant potential equal to 
V = j A. The general solution for ax is 


ax 


Cl 



Q+V 

2Q 


+ C 2 



Q-\/Q^-4uj2 

2Q 


(C.4) 


There are two linearly-independent solutions. When > Q‘^ the first solution describes 
the in-falling wave with the w —?• a; -b ie prescription. When uP < both solutions are 
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Table 1. Geometries and properties of 1, 2, 3, 4-charge black holes in AdS 4 . 



1-charge 

2-charge 

3-charge 

4-charge (RN) 

F(0) 

—6 cosh((?5)/-\/3) 

— 2(cosh ^ -h 2) 

—6cosh(^/-\/3) 

-6 

ZW 
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e'l>/^/3 
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f3/2(f-rQ)l/2 

f(f + Q) 

fl/2(f+Q)3/2 

(r + Q? 
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^ rjj{fH + Q) 

^ rHifH + Q)^ 

{rn + Qf 

^ r{fH + Q)* 

P(r + Q) 

f(f + Q)2 

(f + QP 

rnir + QP 

At 

A rH + Q\ 


, /SClffrr -h 0) ( 1 Q A 

^VQ{rH + Q) 7 rn + Q\ 

pfu -1- Q V f + Q ) 


VAQPh + 1^17 f + QJ 

VTV ^ r + Q ) 

<t> 

#'»{■ +f) 

V-f) 


0 

(7, >5) 


7 = 1, 5 = —1 

1 . 1 

1 = 0 = ^ 

V3 V3 

7 = 0, 5 = 0 

{z,e) 

z = i,e = -i 

z = 1 , 6 ^ 00 

2: 00, —Qjz = 1 

2 00, 9 = 0 


K = Vs, c = -1 

K = 2, C = -1 

K = V3, c = 1 



4 „ 

a = -, /3 = 0 

3 

0 

II 

1 

II 

0 

II 

1 

II 


T 

1 

2V3 

0 

1 

2V3 


T 

(3frr + 2Q)V5hr 

3fH + Q 

3\/rH{rH + Q) 

{rn + Q){3rH - Q) 

AnPru + Q 

Att 

47r 

Airfu 

S 

see^ 

see ^ 

(X T 

const 


Exceptional 

Exceptional ^ 

IR charged, gapless 

IR charged, gapless 


^ They are defined in eq. (2.8). 

The scaling exponent depends on the ensemble (fixing p or fi). See figures 13 and 14. 
For the vh ~ 0 geometry, the conductivity has a gap. See (C.5). 


real and normalizable. As the background extremal solution is singular, this situation seems 
to be in the holographically non-well defined class we have discussed earlier: we need to 
impose an extra boundary condition at the singularity in order to obtain a unique solution. 

However, in this case, analyticity (in w) of the correlator makes the solution unique. 
Indeed, if we analytically continue the > <5^/4 solution to < <5^/4, the solution for 
ttx is unambiguous for all cj. The conductivity is^'^ 


a{u;) = 


i{Q + VQ2 _ 4^2) 


2uj 


(C.5) 


oj-t-u+ie 


assume cu > 0. It Is crucial to shift the pole to the lower half tj-plane: y/—uP' —>■ ^ —[u) ie)^ = 

— it = —iuj. 
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Figure 12. The relation between the free energy and the temperature for 1-charge (7 = -s/S) 
S = —l/y/S), 2-charge (7 = 1, <5 = —1), 3-charge (7 = l/VS, 6 = —1/\/3), and RN black holes in 
AdS4. 




Figure 13. 2-charge black hole. The solid line is for the canonical ensemble (fix p), for which 
T = T/y/p and S = S/p\ the dashed line is for the grand canonical ensemble (fix p), for which 
f = T/p and S = SIp^. 



Figure 14. 1-charge black hole. The solid line is for the canonical ensemble (fix p), for which 
T = T/y^ and S = S'/p; the dashed line is for the grand canonical ensemble (fix p), for which 
f = Tjp and S = S/p2. 
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From (C.5) we observe that the spectrum is gapped and continuous above the gap.^® The 
conductivity has a ^-function at w = 0. The AC conductivity as a function of w is plotted 
in figure 15. We should note that very few exact solutions for the conductivity at arbitrary 
Lj are available for extremal geometries; another example is the 1-charge black hole in AdS^ 
[89], 

1-charge black hole in AdS^. The equation for is 

~CLrp “h 7r~, — 0* 


+ 


r + Q r^{r + Q) 
The effective Schrodinger potential is given by 

16{Q + r) 

The solution with the in-falling wave in the IR is 


Ox — Cl 


r + 




r + 1 


where is the Hankel function. 

From this solution we obtain the conductivity as 


r i -^2 
cr(L 0 ) = i 


( 1 ) 


(2a;) 




(jJ — 


As a; —)■ 0, 


I . 

a = —h (vr -I- 2^7 -I- 2i lna;)a; -|- 


ui 


(C.6) 


(C.7) 


(C.8) 


(C.9) 


(C.IO) 


The spectrum is gapless, and the conductivity has a 5-function at a; = 0. The AC conduc¬ 
tivity as a function of uj is plotted in figure 15. 


D DC conductivity at zero density 


Although we focus on finite density systems, it is illuminating to examine zero density sys¬ 
tems at extremality. Consider a zero density density system whose background is described 
by 


S 


d'^xyj—g 


R--idcpf-V{^) 


(D.l) 


and the dynamics of the gauge field is given by a Z{cl))F‘^ term, where V{(p) and Z((j)) are 
defined by (2.2). At finite temperature, the DC conductivity is finite and can be calculated 
from (5.7) with the IR neutral geometry. The result is 


c^DC = Zh^ r 


K'y 

h 


„-27'5 
'h > 


(D.2) 


where r/i —)• oo is the extremal limit. 

^®This is very reminiscent to the situation studied in [52, 53] in Einstein-Dilaton gravity with a potential 
that asymptotes at strong coupling to the Liouville potential of non-critical string theory and the dilaton 
solution to the well-known linear dilaton solution of string theory. 


44 














2-charge 1-charge 




Figure 15. The left plot is the AC conductivity calculated from the 2-charge black hole in AdS 4 
at extremality, and the right plot is the AC conductivity calculated from the 1-charge black hole in 
AdS 4 at extremality. Both of them have a 5-function at a; = 0. 

• When 7(5 > 0, cjdc —t 0 in the extremal limit, so the system is an insulator. 

• When 7(5 < 0, (Tdc —t oo in the extremal limit, so the system is a conductor. 

For the extremal limit (T —>■ 0) and the DC limit (w —>■ 0) to be compatible the conductivity 
should have a (5-function at a; = 0 when 7(5 > 0, while the conductivity should not have 
a (5-function at cj = 0 when 7(5 < 0. We will examine whether this is true for the gapped 
geometries in the following. 

To obtain the conductivity, we perturb the system (D.l) around the background metric 
(3.1) by 5Ax = aa;(r’)e“*^*, and obtain the equation for ax'- 

j^UZa'x)' + ^ax = ^, (D.3) 

where / = ^JD/B. As noted before, we can ignore the term when we discuss the 
a; —)• 0 limit of the conductivity for extremal geometries at T —)• oo, while we need to do 
the asymptotic match for extremal geometries at T —)• 0. We will discuss the extremal 
geometries at T —?■ oo. In terms of the Schrodinger coordinate ^ in (3.9), the general 
solution of ttx in the IR neutral geometry is 

+ (D.4) 

where uq is defined in (3.18). From (D.3), we observe that there is a radially conserved 
quantity at w = 0: 

n = /Za^, 9^n = 0. (D.5) 

If n —7- 0 as r —)• oo (the IR), the conductivity does not have a (5-function at w = 0. If 
n is constant as r —)• oo, the conductivity has a (5-function at w = 0 and its coefficient 
is proportional to 11(0; = 0). That there is only one normalizable solution for ax requires 
\i'q\ > 1. There are two cases as follows: 
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• Case 1: uq > 1, for which we have <5^ — 2'y6 — 1 > 0. The solution for ax is 

ax ~ The radially conserved quantity 11 evaluated at the IR is 

n = fZa'x —t constant. (D-6) 

• Case 2: uq < —1, for which we have 5^ — 27 ^ — 1 < 0. The solution for ax is 

ax ~ . The radially conserved quantity evaluated in the IR is 

n = jZolx ~ ^ 0 . (D.7) 


We expect that the insulating phase corresponds to the case 1, and the conducting 
phase corresponds to the case 2. We will show that this is true after we excluded the 
holographically unacceptable region |z^o| < 1 have imposed Gubser’s criterion. 

Gapped, conducting phase. The condition that the extremal geometries are at T —)• 00 
limit of the small black-hole branch is > 1. The condition that the system is a conductor 
is 7(5 < 0. These two inequalities imply 36^ — 2'y6 — 3 > 0. Consequently, 


I/O + 1 = 


3(5^ - 27(5 - 3 
2((52 - 1 ) 


> 0 . 


(D. 8 ) 


If we further impose that there is only one normalizable solution, which translates to |i/o| > 
1, we finally obtain i/q > 1. According to the discussion above, this corresponds to the case 
1. The conductivity has a (5-function at w = 0 in the gapped, conducting phase. Note that 
this is not always true without the constraint |i/o| > 1 . 

Gapped, insulating phase. The condition that the extremal geometries are at T —)• 00 
is (5^ > 1. The condition that the system is an insulator is jd > 0. These two inequalities 
give -|- 275 — 1 > 0. Consequently, 


^ (52 + 27(5-1 

2((52 - 1 ) 


< 0 . 


(D.9) 


If we further impose that there is only one normalizable solution, which is |z/o| > 1, we 
obtain z/q < —1. According to the discussion above, this corresponds to the case 2. The 
conductivity does not have a (5-function at w = 0 in the gapped, insulating phase as expected 
from general principles. 

The above conclusion can be obtained in another way. The equations for a^, at w = 0 
can be solved by 


Q-x — Cl C2 


dr' 


(D.IO) 


/o f{T')Z{(^{r')y 

In the IR, the first constant term corresponds to the first term in (D.4), and the second 
term corresponds to the second term in (D.4). In the UV, we have = Ci + C 2 r + • • •, 
and thus the conductivity is 


a = 


C 2 

iuCi 


(D.ll) 


We can see that if only the first solution is normalizable, we have C 2 = 0, and thus there 
is no (5-function at cj = 0; if only the second solution is normalizable, we have C 2 7 ^ 0, and 
thus there is a (5-function at a; = 0. 
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For the finite density systems with axions, we can obtain a solution for Ai similar to 

(D.IO). The hrst term in IT is always dominant in the IR. Equation 5.3 at cj = 0 in the IR 

can be signihcantly simplihed to be (/Z 2 A'^)' = 0, which can be solved as 

~ f dv 

\, = Di + D2j—. (D.12) 

This is the counterpart of (D.IO) for hnite density systems with axions, but the difference 

is that (D.12) is valid only for the IR geometry. 
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